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The density functional theory is extended to account for self-bound systems. To this end the 
Hohenberg-Kohn theorem is formulated for the intrinsic density and a Kohn-Sham like procedure 
for an A'^-body system is derived using the adiabatic approximation to account for the center of 
mass motion. 



PACS numbers: 21.60.Jz,31.15.E- 



o 
o 
o 

> 

o 



> 

On 
O 



O 



X 



Introduction - The formulation of a density functional 
theory (DFT) for self-bound systems is a question of cur- 
rent interest in nuclear physics Nuclear application 
of the DFT through the Skyrme model or using different 
density functionals is the only viable way to study heavy 
nuclei. The basic physical quantity in the DFT is the sin- 
gle particle density. This quantity, being well defined for 
an iV-body system being localized by an external field is 
essentially zero for a freely moving self-bound system [l[ . 
Therefore, when treating finite self-bound many-particle 
systems the DFT should be modified to account for cen- 
ter of mass motion. This modification can be looked for 
along two possible lines, either to add an external po- 
tential, confine the system so that the DFT applies and 
remove the center of mass effects at the end of the calcula- 
tion . The other possibility is to derive the Hohenberg- 



Kohn (HK) theorem 
struct a Kohn-Sham 



density and energy 



for the intrinsic density and con- 
like procedure for calculating the 
In this manuscript we shall follow 
the second path. We shall first derive the HK theorem 
for the intrinsic density. This theorem, as the original 
HK theorem, is a specific case of a more general theo- 
rem by Valiev and Fernando who proved the existence of 
an energy functional for any Hermitian operator [l|, Q . 
Here we rederive this proof using our notation for com- 
pleteness and clarity. Once we have the HK theorem at 
our disposal, we follow the KS procedure, and construct 
a KS-like Schrodinger equation, which can be separated 
into a single particle equation using a Born-Oppenheimer 
type approximation for the "slow" center of mass coor- 
dinate. In this last step one is forced to introduce some 
approximations since the desired KS orbitals break the 
translation invariance. 

For clarity we shall restrict our discussion to a par- 
ticle system interacting via two-body forces which only 
depend on the relative distance between the particles. 
Furthermore, in deriving the KS like equations we shall 
follow the Hartree rather than the Hartree-Fock scheme. 

The HK theorem for the intrinsic density - Consider 
a system of N particles moving under the influence of a 



two-body potential 



JV 



and an intrinsic "one" -body potential 

N 



(1) 



Vin = ^Viniri - R) 



(2) 



where R = is the center of mass coordinate. 

The Hamiltonian describing the internal motion of the 
system is 

H = T„, + U + V,„, (3) 
where the intrinsic kinetic energy operator is given by 



^ „2 p2 

^2 2N 



(4) 



Here P = X^Tpi center of mass momentum. Fol- 

lowing HK we shall assume a nondegenerate ground state 
The intrinsic translational invariant density 



N 



n™(r) = m}_^S{n-R-r)\^) 



(5) 



is clearly a functional of Vin{r). Following the HK 
method we can prove that conversely (up to a constant) 
Vin{r) is a functional of ri,i„(r). We start by assuming 
that another potential d ■„ (r) with a ground state 5*' gives 
rise to the same intrinsic density ni„(r). The two wave 
functions must be different unless Wm — = const, 
since they correspond to two different Hamiltonians H 
and H' — Tin + U + y/„ . Utilizing the variational princi- 
ple one gets 

E' = {^'\H'\^') < {^\H'\^) = {^\H + V^n-V^n\'^f) (6) 



or 
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E' < E 



■ j drn,n{r){v[n{r)-Vin{r)) . (7) 
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Interchanging primed and unprimed quantities, we get 
E<E'+ [ drnM {v^n{r) - v[^{r)) . (8) 



Adding the last two equations we reproduce the famous 
inconsistency 



E + E' <E' + E . 



(9) 



Thus we have shown that Um(r) is a functional of the 
intrinsic density nin{r). Since, Vin(r) fixes H we see 
that the intrinsic iV-body ground state is a functional of 
nin{r). 

Since 5" is a functional of the intrinsic density we can 
define the functional 



FK„(r)] = {■^\T,r, + U\^) 



(10) 



which is a universal functional valid for any number of 
particles as we can repeat the proof for any number of 
particles N which by itself is a functional of nin{r). The 
HK variational principle for the ground state energy 

E[n,n{r)]= dr Ti^n{r)v„i{r) + F[n„,{r)] (11) 



follows trivially from the original proof 

The KS equations - Once it has been established that 
the energy of the finite self-bound many-body system is 
a functional of the intrinsic density the next step is to 
construct extended KS equations for such system. The 
ground state energy of the iV-body system interacting 
via two-body potential u can be written in the form 

^ = \ j j drdr'THn{r)u{r - r')nin{r') + G[?i„i(r)] . 

(12) 

where G is a universal functional of the intrinsic density. 
Following KS we write 



GK„(r)] =Ts[n^„]+E^ 



(13) 



where Ts[nin] is the intrinsic kinetic energy of a sys- 
tem of particles interacting with a "one" -body potential 
Vin[nin], and Exc[nin] is the exchange-correlation func- 
tional. Due to the HK variational principle for the en- 
ergy, Eq. IT^ . we obtain the equation 



dr Sn,n{r) 
subject to the condition 



ST, 



Sni„{r) 



, 



J dr6nin{r) = 



(14) 



(15) 



Here 



and 



SExc[nin] 

5nin{r) 



(17) 



Equations , (fT51) are precisely the same equations one 
obtains starting with the Hamiltonian, 



R) 



(18) 



Unfortunately this Hamiltonian cannot be separated into 
single particle orbitals like the original KS Hamiltonian. 
However, since the center of mass coordinate is a slow 
coordinate we can replace all the center of mass operators 
by their expectation values. In order to do so, let us 
assume that is a smooth function that can be expanded 
in the following manner. 



where 



da'~p{r) = 



dip[r) 
dvadn 



(19) 



(20) 



and summation over the spatial directions a,b — {x, y, z} 
is assumed. Replacing the center of mass terms by their 
expectation values and choosing a coordinate system such 
that (R) = 0, we obtain 



N 



H 



^ 2 



2N 



Y.\^^{T,)^h,RaRb)dl^,^(Ti) 



(21) 



retaining only the leading center of mass corrections. 
The d\a^\vi correction accounts for a linear term in V<p 
which upon summation over all particles is proportional 
to R. Clearly this Hamiltonian is a sum of N single- 
particle Hamiltonians leading to the KS type orbitals in 
the Hartree approximation 



^(r) + i(i?„i?,)a„VW 



= e,?/., . (22) 



To be more concrete let us consider the local density ap- 
proximation (LDA) in which it is assumed that ^Pxcij) = 
(y92,f,(ni„(r)). In this approximation the leading correc- 
tion to the single particle potential due to the center of 
mass is 



1 



(23) 



where 



^{r) = (^,,(r) + y dr'u{r - r')n,„(r') (16) ^iM^) 



dtp.. 



xc a2 



drir 



tpxc 



dnl 



dan{r)dbn{r) 
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dr'u{r - r')dli,n{r') 



(24) 



Here n{r) is the laboratory one-body density given by 

N 

n(r) = 5]|^,(r)r (25) 

i 

Equations HI]) , dTT]) , (I121)-(Il5l) together with 



{RaRb) 



N 



1 ^ 



(26) 



form a closed set of equations to be solved self- 
consistently, just as in the case of a system bounded by 
an external potential. 

The nature of the center of mass expansion depends 
on the system at hand. This point is most easily demon- 
strated in the Hartree picture where the expectation 
value of all the odd moments of Ra can be set to zero 
and the second even moment is given by 



{RaRbRcRd) — 

{{rarb){rcrd) + {rarc){rbrd) + {rard){rbrc)) 



(27) 



From Eqs. (I26p and ((27|) it is clear that the leading orders 
in the expansion pop can be regarded as a Taylor series 
in {r'^)/N, and that in the limit N — > oo the original KS 
equations are recovered. For a collapsing system where 
the size of the system does not depend on the number 
of particles the expansion behaves as a series. 
For non-collapsing systems, such as nuclei, the rms radius 
grows as and the expansion converges much slower 
as 1/ -v^^, unless the r-dependence of the potential scales 
as the rms matter radius, i.e. ip « V'l^/ \^ci). 

Example - As an example for our procedure consider 
a system of N fermions moving in a harmonic oscillator 
potential, 



H 



Using the relation 



N 



N 9 

^ 2 



N 



(28) 



J2ir. -Rf = J2-j-NR' 



(29) 



we can rewrite this Hamiltonian as a sum of a center of 
mass term 



P 1 

Hem. = h -Nuj'^B? 

2N 2 



and an internal Hamiltonian written in the form (1181) 



N 2 p2 1 ^ 

^'^"-Ef -^ + ^E->^-«)'- (31) 

i ___ i 

Following Eq. (PT|) . we evaluate the center of mass cor- 
rections 



\{RaRb)dlMr.) = l{R')c^' , (32) 



and get the approximated internal Hamiltonian 



^» ^ E Y - ^ - l^-'(R') + E l-'r^^ • (33) 



(30) 



Comparing (|33|) with ([30)) it is clear that for the harmonic 
oscillator case the corrections to the internal Hamiltonian 
ensure the right cancellation of the center of mass energy. 

Discussion - It is already known for some time that 
the HK theorem can be generalized to any Hermitian 
operator, including the intrinsic density. The challenge 
however is to find a useful KS like procedure that re- 
duces the self-bound many-body system into a set of 
single quasi-particle orbitals. From the onset it is clear 
that such orbitals break the translational symmetry, and 
therefore some approximations are called for. In this note 
we have demonstrated that starting with the HK theo- 
rem for the intrinsic density one can reach this goal by 
treating the center of mass coordinate as an adiabatic 
variable. This approximation recovers the KS equation 
in the limit N — > oo and includes center of mass recoil 
effects in the orbitals equations and to the ground state 
energy. We have demonstrated that for the case of A^- 
particles moving in an harmonic-oscillator potential this 
procedure yields the right center of mass correction. 

We argue that the importance of the center of mass 
corrections might vary for different systems. Regardless 
whether the effect is small or large the current approach 
helps bridge the gap between the DFT and its application 
to finite self-bound systems. 
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